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Abstract 

This paper is devoted to the characterization of differentially flat nonlinear systems. Implicit 
representations of nonlinear systems, where the input variables are eliminated, are studied in the 
differential geometric framework of jets of infinite order. In this context, flatness may be seen 
as a generalization of the property of uniformization of Hilbert's 22nd problem. The notion of 
Lie-Backhmd equivalence, introduced in [16 , is adapted to this implicit setting. Lie-Backlund 
isomorphisms associated to a flat system, called trivializations, can be locally characterized in 
terms of polynomial matrices of the indeterminate 4k ■ Such polynomial matrices are useful to 
compute the ideal of differential forms generated by the differentials of all possible trivializations. 
We introduce the notion of a strongly closed ideal of differential forms, and prove that flatness 
is equivalent to the strong closedness of the latter ideal. Various examples and consequences 
are presented. 

Keywords. Nonlinear system, implicit system, manifold of jets of infinite order, Hilbert's 22nd 
problem, polynomial matrices, ideals, differential forms, differential flatness, fiat output. 

Introduction 

Differential flatness, or more shortly, flatness, is a system property introduced more than ten years 
ago in [261 [T3]. 

Let us briefly state an informal definition, that will be made more precise later: given a nonlinear 
system 

x = f(x,u) (1) 

where x = (x\, . . . , x n ) is the state (belonging to a given smooth n-dimensional manifold) and 
u = (u%, . . . , u m ) the control vector, m < n, the system ([TJ is said to be locally (differentially) flat 
if and only if there exists a vector y = (y\, . . . , y m ) such that 
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• y and its successive time derivatives y,y, . . . , are locally independent, 

• y is locally a function of x, u and a finite number of time derivatives of the components of u, 

• x and u can be locally expressed as functions of the components of y and a finite number 
of their derivatives: x = <po(y,y, . . . , y^), u = ipi(y,y, . . . ,y^ a+1 '), for some multi-integer 
a = (ai, . . . , a m ), and with the notation = ( d dt L^ 1 , • • • , d J^^T ) • 

The vector y is called a _/?ai output. 

This concept has inspired an important literature. See [131 [27] for surveys on flatness and its 
applications. Let us just mention that flatness provides significant simplifications for the motion 
planning problem or for several aspects of feedback design. 

Various formalisms have been introduced: finite dimensional differential geometric approaches 
[UJ H21 [43] , differential algebra and related approaches [HI EJ [21] , infinite dimensional differential 
geometry of jets and prolongations [HI [29J, [3U |3T| [38]. 

Note that flatness may be traced back to Hilbert and Cartan |19]l5j. In fact, using the definition 
presented in this paper, with, in place of ([I]), the set of n — m implicit equations ([2]) where the 
control variables u are eliminated, flatness may be seen as a generalization in the framework of 
manifolds of jets of infinite order of the uniformization of analytic functions of Hilbert 's 22nd 
problem [18J , solved by Poincare [32J in 1907 (see [3 J for a modern presentation of this subject and 
recent extensions and results). This problem consists, roughly speaking, given a set of complex 
polynomial equations in one complex variable, in finding an open dense subset D of the complex 
plane C and a holomorphic function s from D to C such that s is surjective and s(p) identically 
satisfies the given equations for all values of the "parameter" p G D. In our setting, C is replaced 
by a (real) manifold of jets of infinite order, a flat output yi, . . . , y m plays the role of the parameter 
p and s is the associated Lie-Backlund isomorphism s = (<£o, fi, <fi, ■ ■ •) with ifo and ipi defined 
above. 

In the framework of linear finite or infinite dimensional systems, the notions of flatness and 
parametrization coincide as remarked by [36J [37], and in the behavioral approach of [33J, flat 
outputs correspond to latent variables of observable image representations |44] (see also [12] for a 
module theoretic interpretation of the behavioral approach). 

The characterization of differentially flat systems has aroused many contributions [21 El [TT1 
|2T| [28| [30l [35l [38| 1391 W2\ 03]. Though general necessary and sufficient conditions exist (see e.g. 
[21 El [30]), they don't provide a fully computable set of conditions. More precisely, [2] gives an 
algorithm to compute a basis of the cotangent module, called infinitesimal Brunovsky form, and 
further integr ability conditions are needed to deduce flat outputs. Recently, [9] has improved this 
result but still requiring the use of deformations, thus assuming that the system depends on a small 
parameter. In [30] the author proposes to express all the differential relations between the system 
variables x, u and a candidate flat output y and use Cartan-Kahler theory. 

We adopt here the formalism of manifolds of jets of infinite order [Tl [TB ] [24" 1 134 1 145] and, as previ- 
ously mentioned, in place of systems in explicit form (pQ), we consider (equivalent) implicit systems 
obtained from ([1]) by eliminating the input vector u. This representation has the major advantage 
to be naturally invariant by endogeneous dynamic feedback (see |16] ) . We adapt the notions of Lie- 
Backlund equivalence and Lie-Backlund isomorphism in this context and show, after restricting to 
the category of meromorphic functions, that the flatness property is naturally described in terms 
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of polynomial matrices and differential forms deduced from the variational system equations. For 
a detailed presentation of polynomial rings and non commutative algebra, the reader may refer to 
[1U\ 123] and for exterior differential systems to [I] . 

Though our results show some parallelism with those of [21 [9] , in particular concerning the study 
of variational properties, they differ from the latter by the fact that, as previously announced, they 
are invariant by endogeneous dynamic feedback. In addition, they exploit different ideas that may 
be seen as an extension of the linear approach presented in [25] to nonlinear systems and they 
provide effective flatness conditions as attested by the examples of the last section. 

The paper is organized as follows: the first Section is devoted to the basic description of implicit 
control systems on manifolds of jets of infinite order. In Section 2, we extend the notions of Lie- 
Backlund equivalence and Lie-Backlund isomorphism to the implicit system framework. Section 
3 deals with the presentation of some variational properties of flat systems. The necessary and 
sufficient conditions for flatness are stated in Theorems [2] and [3J of Section 4 and some consequences 
are presented. Section 5 is then devoted to examples and some concluding remarks are given. An 
appendix on polynomial matrices and their Smith decomposition is provided. 



1 Implicit control systems on manifolds of jets of infinite order 

Given an infinitely differentiable manifold X of dimension n, we denote its tangent space at an 
arbitrary point x G X by T X X, and its tangent bundle by TX = U^ex T X X (identified with the 

vector bundle TX ^ X). Let F belong to C°°(TX; M n ~ m ), the set of C°° mappings from TX to 
W a ~ m . In the sequel, we will call smooth any function of class C°° in all its variables. 
We consider an underdetermined implicit system of the form 

F(x,x) = (2) 

regular in the sense that rank = n — m in a suitable open subset of TX. 

Remark 1 Note that any explicit system of the form x = f(x,u) with f smooth, f(x,u) G T X X 
for every x G X and u in an open subset U of W l , and rank (M^j = m in a suitable open 
subset of X x U, can be locally transformed into H^); permuting the lines of f , if necessary, such 
that the m last lines of f are locally independent functions of u, and still noting x the permuted 
vector (this abuse of notations being unambiguous) , thanks to the implicit function Theorem one 
gets u = fi(x,x n - m+ i, . . .,x n ), /x smooth, and ii = fi(x, fj,(x, x n - m +i, ■ ■ ■ ,x n )) for i = 1, . . . ,n - 
m. Thus, setting Fi(x,x) = ±i — fi(x, fx(x, x n - m +i, . . . , x n )), i = 1, . . . , n — m, the system is 
in the implicit form and since — diag^^—^, G^\, G being the matrix made of the entries 
SfcLi J^~§!f f or ^j = n—m+1, . . . ,n, which don't depend on (xi, . . . ,x n - m ), we have rank (§f-) = 
n — m. 

Conversely, (0|) can be transformed into the explicit system x = f(x,u) with f(x,u) E T X X 
for every x G X a and every u in an open subset U of M. m with rank = m: the rank con- 

dition rank (^f-) = n — m and the implicit function Theorem yield Xi = fi(x, x n - m +x, . . . , x n ), 
i = 1, . . . , n — m, with fi smooth for i = 1, . . . n — m, and, setting x n _ m+ j = Uj for j = 1, . . . , m 
(thus u = (ui, . . . , u m ) G W m ), we finally get x = f(x, u) with / n _ m+ j(x, u) = Uj for j = 1, . . . , m, 
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and, by definition, f is a smooth vector field on X for every u in an open subset U ofM. m , with 
^ = I m the identity matrix of M. m . 

Clearly, the underdetermined character of (|2|) is expressed by the rank condition rank(^£) = 
n — m, which means that the system effectively depends on m independent control variables. 

A smooth vector field f that depends, for every x € X, on m independent variables u € M m 
with rank ($£) = m in a suitable open set of X x R» and that satishes F(x, ,(*, „)) _ for every 

u G U, is called compatible with (J2]). 

We also introduce the following definition: 

Definition 1 We say that system f2]j admits the equivalent representation (TJj) at (xo,uo), if and 
only if every smooth integral curve of fip is a smooth integral curve of (flj) passing through (xq,xq), 
with x = f(x ,u ). 

System {TP admits the locally equivalent representation Hj) if and only if |7]) is equivalent to (fj|) at 
every point (xq,uq) of an open subset of X x W 71 . 

In [16] (see also [34J where a similar approach has been developed independently), infinite 
systems of coordinates (x,u) = (x, u, u, . . .) have been introduced to deal with prolonged vector 
fields J(x,u) = Y% = i fi( x ^ u )wi + Y!j=i Efc>o n f +1) ^T' for ex P licit systems i = f(x,u). 

Here, we adopt an external description^ of the prolonged manifold containing the solutions of 
P|): we consider the infinite dimensional manifold X = f X x = f X x W 1 x K. n x . . . made of an 
infinite (but countable) number of copies of M. n . To endow X with a suitable topology, we define 
the continuity and differentiability of functions from X to R as follows: 

Definition 2 We say that a function cp from X to ~K is continuous (resp. differentiablej if cp 
depends only on a finite (but otherwise arbitrary) number of variables and is continuous (resp. 
differentiable) with respect to these variables. 

Thus X is endowed with the coarsest topology that makes projections to products of X and a finite 
number, k, of copies of W 1 , namely X x K fcn , continuous for any k £ N. This topology can be 
identified with the product topology of the infinite product X x W n x W 1 x . . . 

C°° or analytic or meromorphic functions from X to M are then defined as in the usual finite 
dimensional case since they only depend on a finite number of variables. 

We assume that we are given the infinite set of global coordinates of X: 

— V*<1} • • • j J "ni • • • ! ^li • • • i ^nj • ■ ■ ) -^n > • • •/ v°7 

and we endow X with the so-called trivial Cartan vector field [2U [45] 

^EE^ 1 '^. (4) 

i=l j>0 ox i 

1 in the sense that the system manifold, namely the set of points x — (x, x,x, . . .) such that F(x, x) — 0, is described 
by means of the larger manifold X = X x RJ£,. 
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so that x\ k ^ = 4^ for every i = l,...,n and k > 1, with the convention ir^ = Xj, where 
A x (i) 4e. f L^Xi = = x^ , and where 



»=i i>o ax i 

is the Lie derivative of a function (/? € C°° (X, R) along the vector field T£ (this series having only 
a finite number of non zero terms according to Definition [2]) . Therefore the Cartan vector field 
acts on coordinates as a shift to the right and X is called manifold of jets of infinite order (see 

EE23H5I). 

From now on, x y, ... stand for the sequence of jets of infinite order of x, y,. . . 
The implicit representation ^) , as opposed to (pTJ) , is invariant by endogeneous dynamic feedback 
(see [HI Section IV] for a precise definition): 

Proposition 1 Consider the explicit system |7]J with the endogeneous dynamic feedback 

u = a(x,z,v), z = b(x,z,v). (5) 

If the system FT]) admits the locally equivalent implicit representation then the closed-loop 
system 

x = f(x, a(x, z, v )), z = b(x,z,v) (6) 

also admits the locally equivalent implicit representation In other words, system is invariant 
by endogeneous dynamic feedback. 

Proof. We assume that ([I]) admits the locally equivalent implicit representation ([2]) and that 
the endogeneous dynamic feedback (|5|) is given, with z remaining in Z, a given finite dimensional 
smooth manifold. By definition of an endogeneous dynamic feedback, the closed-loop system ([6]) is 
Lie-Backlund equivalent to ([T]). Thus, there exist finite integers a and f3 and locally onto smooth 
mappings and such that 

(x,z,v) = $(x,u, . . . (x,u) = W(x,z,v, . . . (7) 

According to Remark [TJ there exists a locally defined smooth mapping fx such that (H|) is locally 
equivalent to ([2]) with u = fx(x,x). Denoting by = // fc )(x,x, . . . ,x^ k+1 ^) and $(x, . . . = 
$>(x, [i(x,x), . . . ,fjS a >{x, . . . ,x^ a+1 ^)), we immediately deduce that for every smooth local integral 
curve t i— > x(t) of ([2]) passing through an arbitrary point (xo, xq), using the first relation of ((7J), t t— » 
(x(t),z(t)), given by (x(t), z(t),v(t)) = $(x(£), . . . ,x^ a+1 \t)), is also a smooth local integral curve 
of the closed-loop system © passing through (xo, zq) such that (xq, zo,vq) = $(xq, • • • , Xg a+1 ^). 

Conversely, if 1 1— > (x(i), z(t),v(t)) is a smooth local integral curve of ([6]) passing through (xo, zq), 
by the second relation of 0, we get that (x(t),u(t)) = *Zf(x(t), z(t), v (t), . . . ,v^(t)) is a smooth 
local integral curve of (PQ) and, according to the assumption, t \— > x(t) is also a local smooth integral 
curve of ([2]) passing through the corresponding point (xq,xq), which achieves the proof. ■ 
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Definition 3 A regular implicit control system is defined as a triple (X, t%, F) where 
X = X x R^, t x is the trivial Cartan field on X, and F G C°°(TX; M n_m ) satisfies rank (§§) = 
n — m in a suitable open dense subset ofTX. 

Note that if t i— > (x(t),x(t)) is an integral curve of P|), then i i— > is an integral curve of 
L^ X F = for every k. Therefore, every integral curve of the system (X, tx, F) lies in the set of x 
such that L^F = for every k. 

The system (M™,r m ,0), where we have noted r m the trivial Cartan field of M™, is called the 
trivial system of dimension m since it corresponds to the void (unconstrained) implicit system 
= 0. Note that since n — m = 0, the rank condition on 4j 3 which in this case is the matrix whose 
entries are all identically 0, is globally satisfied. The absence of equations relating the coordinates 
of makes this system also an explicit one and, in the notations of [16], it corresponds to the 
trivial explicit system (M.™,T m ), which justifies its name. 

2 Lie-Backlund equivalence for implicit systems 

Let us slightly adapt the notion of Lie-Backlund equivalence^! of [16] in our implicit control system 
context: 

Let us consider two regular implicit control systems (X, t x , F), with X = X x RJ^, dim X = n, 
Tx the associated trivial Cartan field and rank = n — m, and (2),r*g,G), with 3)=7x ffi.^, 

diml" = p, T<g the associated trivial Cartan field and rank (jjijf) = P ~ 1- 

Set X = {x G 3L\L k Tx F{x) = 0, VA; > 0} and 2) = {y 6 Z)\L^G(y) = 0, VA: > 0}. They are 
endowed with the topologies and differentiable structures induced by X and 2) respectively. 

Definition 4 We say that two regular implicit control systems (X,Tx,F) and (fZ),Tsn,G) are Lie- 
Backlund equivalent (or shortly L-B equivalent) at (xQ,y ) G Xo x 2)o if an d on ^V if'- 

(i) there exist neighborhoods Xo and °/ #o in Xo and y * n ?)o respectively and a one-to- 
one mapping <I> = (ipo, ipi, . . .) G C°°(yo; ^o) satisfying Q(y ) = xq and such that the trivial 
Cartan fields are & -related, namely <J>*Tsrj = tx; 

(ii) there exists a one-to-one mapping ^ = (ipo,ipi, . . .) G C°°(Xo;yo) such that ^(xq) = y and 

The mappings <3? and ^ are called mutually inverse Lie-Backlund isomorphisms at (xo,y ). 

The two systems (X,t%,F) and (%},Ttg,G) are called locally L-B equivalent if they are L-B 
equivalent at every pair (x,9(x)) = ($>(y),y) of an open dense subset Z of Xo x 2)o> with $ and ^ 
mutually inverse Lie-Backlund isomorphisms on Z. 

2 In the terminology of [16], different names have been introduced, which, in the author's opinion, are poorly 
matched: differential equivalence (corresponding to endogeneous transformations and ^-related Cartan fields) and 
Lie-Backlund equivalence (including time scalings into the previous endogeneous transformations by replacing Cartan 
fields by Cartan distributions). In order to stress that they are two faces of the same coin, we propose to use Lie- 
Backlund equivalence (resp. orbital Lie-Backlund equivalence) in place of differential equivalence (resp. Lie-Backlund 
equivalence) and Lie-Backlund isomorphisms (resp. orbital Lie-Backlund isomorphisms) in place of endogeneous 
transformations (resp. Lie-Backlund isomorphisms). 
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It is easily seen that two systems x = f(x, u) and y = g(y, v) with / compatible with (X, t%, F) 
and g compatible with (%),Tsg,G), are differentially equivalent in the sense of [16] (or L-B equiv- 
alent as proposed in footnote 2 ) at the pair ((xq, uq, uq, . . .), (yo, vq, vq, . . .)), with uq such that 

= f(xo,Uo) and uo such that yo = g(yo,vo), if and only if (X,t%,F) and (2),T2),G) are L-B 
equivalent, in the sense of Definition [H at (x~o,y ) with x = (xq,xq, . . .) and y = (yo,yo> ■ ■ • ), 
using the construction of Remark [TJ let <3? and satisfy (i) and (iij. For every g compati- 
ble with (2),Tgj, G) and for y £ y we have G(y, g(y,v)) = for all v in a suitable open sub- 
set of W 1 . Since, by assumption, x = 3>(y) G Xo C Xo, with $ = ((pg, (fx, . . .), we have 

x = tpo(y, g(y, v), ^f(y, v, v), . . .) = (/5o(y,v), and x satisfies L^ x F(x) = for all A; > 0. Let 7>o 
be such that y = g(y ,v ). Thus ^ = ®(yo, g(yo, v ), -^g(yo,v ,v ), ■ . •)• For every / compatible 
with (X,r x ,F) one has x = f{x,u) for u = fi{x,x), or u = ^{(p a {y,v), ^<p {y,v)) = (pi{y,v). Using 
$*r«g = tx yields 

f ( x > u ^\(w(v,v),Mv,v)) = Lr x x \(My,v)m(y,v)) = L ^jMy, LgV, I%V, • • •) 
j>oi=idyl v j>ai,k=i H ° Vk j,i>ok=n=i dVi dvh 

k>0 yK l>0 k=l ov k 

Analogously, u = ju = L T ,^(x) = |^(y,«) = £? =1 ft (v,«)^ + E,>o E?=i ^If) = 

(<Pi)*9, which proves that (x,u) = $(y,v), with $ = (<po, <pi, . . .), and / = <I>*y, defining uo 
by (xo,tio) = ^(yoi^o)- Symmetrically, we have (y,v) = ^(x,u) with \& = (tpo,tpi, . . .) and 
■4>q{x,u) = ip (x, f(x,u), %(x,u, u), . . .), = v{i>o(x, u), £i> (x,u)), v = u(x,x), and thus 

$'s inverse is $ with (yoj^o) = ^(^o^o), and g = which proves that the explicit systems 

(X x R™, /) and (Y x R|o,y), for every Cartan fields / compatible with F = and g compatible 
with G = 0, are locally L-B equivalent at {xq,Uq), (yo,vo) for uo and uo suitably chosen, their 
choice depending on / and g. The proof of the converse follows the same lines. 

An easy consequence of this definition is that L-B equivalence preserves equilibrium points, 
namely points y = (t, y, 0, 0, . . .) (resp. x = (t, x, 0, 0, . . .)) such that G(y, 0) = (resp. F(x, 0) = 0). 

The following result is easily adapted from |16j : 

Proposition 2 If two regular implicit control systems (3C,tx,F) and (%),Tsg,G) are locally L-B 
equivalent, they have the same coranks, namely m = q. 

3 Flatness and variational properties 

First recall from [16] that a system in explicit form is flat if and only if it is L-B equivalent to 
a trivial system. The reader may easily check that this definition is just a concise and precise 
restatement of the definition given in the Introduction. The adaptation of this definition in our 
context is obvious: 
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Definition 5 The implicit system (%,tx,F) is flat at (xo,y ) € Xo x if and only if it is L-B 
equivalent, at (xo,y ) £ Xo x M™, to the m- dimensional trivial implicit system (M™,T m ,0). In 
this case, the mutually inverse L-B isomorphisms $ and are called inverse trivializations, (or 
uniformizations, in reference to Hilbert's 22nd problem). 

Otherwise stated, ([2]) is fiat at (xo,y ) if ^he local integral curves t i— ► x[f) of 
^ around xq are images by a smooth one-to-one mapping $, satisfying x = &(y ), 
of arbitrary curves in the coordinates y = (yi, . . . , y m , y\, . . . , y m , . . . , , . . . , y$ ,•••,) 
around y . In other words, for every curve t i— > y(i) in a suitable time interval J, 



z(t) = (a:(t), x(t), x(t), . . .) = $(y(t)) = {y (y{t)), y>i(y(i)), p 2 (!/(*)), • • •) belongs to X for all t G J 



and thus F(x(t), x(t)) = 0. Conversely, if t \— > x(i) is an integral curve of F(x,x) = 0, there exists 
a curve t i-> y(i) in C°°(J;]R m ) such that y(t) = (y(t),y(t), . . .) = V(x(t)) = (il) (x(t)), ij)i(x(t)), . . .) 
for all i G J, namely L Ty y(i) = L TX ^/(x(t)) for all t G 1 Recall that $ (resp. \P) depend only on a 
finite number of derivatives of y (resp. x). 

The extension of this definition to local flatness is straightforward. 

Trivializations may be characterized in terms of the differential of F. A basis of the tan- 
gent space T^X of X at a point x € X consisting of the set of vectors { ti s \i = 1, ... ,n, j > 0}, 

dxy> 

a basis of the cotangent space Ti-X at x is therefore given by {dx[^\i = l,...,n,j > 0} with 
< dxf j — ^jj >= Si^Sj^i, $i,k being the Kronecker symbol (i.e. 5ij = 1 if i = j and 8^ = other- 



Since smooth functions depend on a finite number of variables, their differential contain only a 
finite number of non zero terms. Accordingly, we define a 1-form on X, an open dense subset of X, 
as a finite linear combination of the dxf \ with coefficients in C°°(X;1R), or equivalently as a local 
C°° section of T*(X x R n ) (see e.g. [35]). The set of 1-forms on X is denoted by A*(X). Clearly, 
the dFi's are elements of A 1 (X). 

Note that the shift property of = L Tx on coordinates extends to differentials by defining 
^dxi as: ^dxi = dii = d-^Xi. More generally, we define the Lie-derivative of a 1-form co along a 
vector field v on X, which is a 1-form on X, denoted by LyLd, as in the finite dimensional case by 
the Leibnitz rule: 



for every vector field w on X, where \v, w] is the Lie-bracket of v and w. Clearly, if u = dxi, v = tx 
and w = ^7 , we recover the previous formula. 

If $ is a smooth mapping from 2} to X, the definition of the image by $ of a 1-form is the 
same as in the finite dimensional context: if w G A 1 (X), uj(x) = Yl i>o Ylt=l \ its 



wise). The differential of Fi is thus given by 




{Ljjuj, w) = Ljj (u, w) - (uj, [v, w]) 



finite 



(backward) image is the 1-form on 2) defined by 




(') 



(9) 
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where <pij is the i,j-th component of <£, namely = (pij(y). 

Note again that, since the functions (fij depend on a finite number of variables, the 1-form $*cu 
contains only a finite number of non zero terms and, according to x^p = tI^, we have 

- E a -m *P = = *4<* - ^ ( E ^4" 



kl u yk \ kl 




(10) 



Theorem 1 TTie system {X x R^,, tx,F) is flat at (xo,y ) with xq G Xo and y € R™ i/ and on/?/ 
if there exists a locally smooth invertible mapping $ from R™ to Xo, with smooth inverse, satisfying 
$(y~o) = xq, and such that 

$*dF = 0. (11) 

Proof. Necessity: If the system (X x R^, tx,F) is flat at (x , y ), we have F(ipo(y), (fi(y)) = for 
all y in a neighborhood of y in RJJ. For every A in a given interval [0, Ao[ of R and a sufficiently small 
time interval J containing 0, consider a smooth trajectory 1 1— > y A (i) i n a bounded neighborhood of 
yo in R m , such that sup{||y A "^(i)|||j > 0, t G 3, A G [0, Ao[} is finite, where || • || denotes the Euclidian 
norm of R m , and set ^^-(t)u_ = C(^)> which exists on J by Ascoli's Theorem. Next, we set 

x x (t) = *(y A (i)) for t G 3 and A G [0, A [. We indeed have F(<p (y x (t)),(pi(y x (t))) = for all 
t and A in their respective intervals and thus, differentiating with respect to A at A = 0, we 
§ et m T§Ki7o(t))C(t) + M^r(ya(t))((t) = in J. At time t = 0, noting y A (0) = y A>0 , we have 
S^(%,o)C(0) + H^(yo,o)C(0) = 0, and this expression is valid for every y 00 in a neighborhood 

of y and C(0) G T %0 R^. We have thus proved that the 1-form ff^dy + §f ^dy = $*dF 
vanishes on Tj^ MJJ for every y 00 in a neighborhood of yo, and therefore is identically zero in a 
neighborhood of yo, which proves (jlip . 

Let us prove the sufficiency. Assuming that there exists a locally smooth invertible mapping 
$ = (fo,fi, ■ ■ •) € C°°(R™;X) satisfying ([TTJ with $(y ) = #o, the 1-forms <I>*dFj are obviously 
closed since they are the differentials of the functions F o i = 1, . . . , n — m. Thus (jlip implies 
that Fi((po(y), <£>i(y)) = Q, i = 1, ...,n — m, with Cj arbitrary constants. But since xo G Xo 
and x = $(y ), we have F(x ,»o) = and Cj = F(<£> (yo); ¥>i(l7o)) = ^i(^o,»o) = 0, for i = 
1, . . . , n — m. Then, setting x = 3>(y) = (<£>o(y)> ¥>i(y)> • • •)> we get that x = y?o(y) (which depends 
on y and a finite number of its derivatives) satisfies F{(x,x) = 0, i = 1, . . . ,n — m and that 
x = L Tx x = L Ta) (y9o(y) = ¥>i(y)- Following the same lines for all derivatives of x, we have proved 
that <1?*T2) = t%. Finally, since is invertible with C°° inverse \P = (^Q,tpi, . . .), it is immediately 
seen that y = ^(x) (and therefore y = tpo(x) only depends on a finite number of derivatives of x) 
and that ^*t% = rm, which proves the sufficiency and the proof is complete. ■ 
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4 Flatness necessary and sufficient conditions 



according to (jlOp . we have: 



We now analyze condition (jlip in more details with the (mild) restriction that F is meromorphic 
on TX. This restriction is motivated by the use of algebraic properties of polynomial matrices and 
of modules over a principal ideal ring of polynomials, this ring being itself formed over the field of 
meromorphic functions, as will be made clear immediately. 

We also restrict the inverse trivializations <3? and of definition [5] to the class of meromorphic 
functions. 

In matrix notations and using indifferently ^ for L Tx or L T ^ (the context being unambiguous), 



dF difp dF dipi \ ^ 
dx dy U) 8x QyU) J Vi 

9F 9<Po dv® + —-( 9<P0 dv® 

dXQyO) ^ ^ didt (?) ^ 

8F_ dF_d_\ (sTX^ d( PQ dj d 

Introducing the following polynomial matrices where the indeterminate is the differential oper- 
ator |: 




with P(F) (resp. P(<£>o)) of size (n — m) x n (resp. n x m), (TTT1) reads: 



or equivalently 



$*dP|_ = P(F)| S(5) P( V0 )| 5 <fo = 0. (13) 

* m w\ m = p W\* p to)\ m d y = °- ( 14 ) 

Clearly, the entries of these matrices are polynomials of the differential operator 4z whose coefficients 
are meromorphic functions from X to R. We denote by R the field of meromorphic functions from 
X to R and by the principal ideal ring of polynomials of 4 = L T% with coefficients in K. 

We may also consider the field of meromorphic functions from 2) = R™ to R. In this case, the 
notations and with ^ = L T?) , will replace the previous ones. 

Recall that &[4z] is non commutative, even if n = 1, as shown by the following example: 
denoting, with abusive notations, by x the 0th order operator a ^ xa for every o 6 ^, we have, for 
a ^ 0, (4r ■ x — x ■ 4f) (a) = xa + xa — xa = xa ^ 0, or ^ ■ x — x ■ = x. 

For arbitrary integers p and q, let us denote by M Pi9 [^] the module of p x q matrices over J^[4] 
(see e.g. |10} 123] for a detailed presentation of modules over non commutative rings). Recall that, 
since in general the inverse of a polynomial is not a polynomial, for arbitrary p£N, the inverse of 
a square invertible matrix of doesn't generally belong to Matrices whose inverse 
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belong to M PiP [^] are called unimodular matrices and their set is denoted by U p [j|]. It forms a 
normajl subgroup of the group of invertible matrices of M P)P [^]. 
Here, P(F) G M 

n—m,n[df] and, according to the Annex [A] (see also [TO]) it admits a Smith 
decomposition (or diagonal reduction), given by 

VP{F)U = (A,0 n _ m , m ) (15) 

with n _ mim the (n — m) x m matrix whose entries are all zeros, V G U n _ m [^], U G U- n [^] and 
A G M n _ mj „_ m [^] a diagonal matrix whose entries dij divide djj for all < i < j < n — m. 
Moreover, A is unique. 

Consider a point x G X and its projection x on the original manifold X. Let us denote by 
£ = (£i> ■ ■ ■ > £n) a basis of the tangent space T X X. We also denote by [£] the ^[^]-module generated 
by £i by [P(F)$] the submodule generated by the lines of P(F)£ and by M. the quotient module 
[£]/[P(F)£]. We call M the variational module of ((2J) at x. 

Clearly, for the trivial system 2) = , its variational module M.<x) at an arbitrary point y G 
reduces to the free ^[^J-module [C] generated by ( = (Ci, • • • , Cm), a basis of the tangent space 
T y M. m at y. 

According to |11| . A4 can be uniquely decomposed into the following direct 
sum: Ai = T © J-, where T is torsion and T is free. It is immediate to ver- 
ify that T = {£ G M\5C = (Ci> • • • > Cn-m) s.t. £ = ( °" l ' n_m ) C} = im 7r n _ m where 
T^n—m i s the projection operator 7r n _ m 

(£i, •••,£«) = (0) • • • ,0,£ m+ i, . . . ,£ n ) and 

T = {£G.M|(I m ,0 

m,n— m) £ 7^ and A£ — 0} — ker A D ker 7r n _ m . Therefore we have 
dim^ .M = n — m + dim^ (ker A n ker 7r n _ m ) . 
We define the rank of P(F) by rank (P(F)) = dim^ A4 = n — m + dim^ (ker A n ker 7r n _ m ). 
Definition 6 Given a matrix M G 3Yt Pj? [-§|], we say that M is hyper-regular if and only if its 

Or 



Smith decomposition leads to (I p ,0 P) g_ p ) if p < q, to I p if p = q and to ( 9 ) if p > q. 



Clearly, the rank of an hyper-regular matrix is equal to min{p, q}. Moreover, a square matrix 
M G Mp iP [-37] is hyper-regular if and only if it is unimodular: M G 1X P [^]. 

4.1 Flatness and controllability 

Recall from [11] that the variational module of ([2]) at x is said to be controllable if and only if it is 
free. 

Proposition 3 If system (0|) is locally flat at xq, its variational module at every point of a neigh- 
borhood ofxQ is controllable and P(F) is hyper-regular in this neighborhood. 



3 A- 1 U p [4 i ]A = U p [^] for every invertible A £ M p , p [£ t 
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Proof. Let (x,y) belong to a neighborhood of (xo,y ) in Xq x with x = ${y), <& being a 
trivialization. We denote as before ipo the first component of <3?, namely x = fo(y), and we consider 
the variational module A4 = [£]/[P(P)£] at x. Since $ is a trivialization, it is not difficult to check 
that the Jacobian matrix ^pr has rank n in the considered neighborhood and that for every element 

C G A^m ; there exists £ = P((/?o)C G the converse being also true since the inverse trivialization 
^ enjoys the same properties. Moreover, since L%. dtpo = d(L^yj ) = difk for every k, it results 
that Ai is isomorphic to M<g which is free. Hence A4 is controllable. 

The variational module A4 at x being controllable, we have T = ker Anker 7r n _ m = {0}, A being 
the diagonal matrix resulting from the Smith decomposition of P(F). But if £ G ker A n ker 7r n _ m , 
according to the fact that the diagonal elements d^i of A are polynomials of ^, £ satisfies the set 
of differential equations di^ = 0, i = 1, . . . , n — m. And since d\ \ divides da for i = 1, . . . , n — m, 
di,i£ = implies £ = if and only if d\^ G R. Assuming that d^, G .ft and applying the same 
argument to dj+i^+i, we obviously get that di t i G .ft for all i = 1, . . . ,n — m which proves that 
all the diagonal elements of A belong to .ft. It is thus straightforward to modify the unimodular 
matrices U and V such that d^j = 1 for all i, which proves the Proposition. ■ 



4.2 Algebraic characterization of the differential of a trivialization 

From now on, we assume that P(F) is hyper-regular in a neighborhood Xo of xo G 3Co- 
In other words, there exist V and U such that 

VP{F)U= (/ m ,0 n _ m>m ). (16) 

C7 and V satisfying (|16p are indeed non unique. We say that U G R — Smith (P(F)) and 
F G L - Smith (P(P)) if they are such that VP(F)U = (J m ,0). 
Accordingly, if M G M 

n,m[jr]) is hyper-regular with rank(iV/) — m < n, we say 
that VeL-Smith(M) and W G R — Smith (M) if V G U n [|] and W G U m [|] satisfy 

VMW =( T q 

In this framework, the set of all polynomial matrices Pifo)^^ £ '^■n,m{^\ satisfying (fT3l) . or 
(I14p . can be completely characterized. We first solve the matrix equation: 

P(F)Qdy = (17) 

where the entries of G G M njm [^] are not supposed to be gradients of some function ipQ-. 

Lemma 1 The set of hyper-regular matrices G M„ )m [j] satisfying fliTj ) is nonempty and given 
by 

e = u( \ m ' m ) w (is) 



with U G R - Smith (P(P)) and W G U m [4] arbitrary 
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Proof. Using Theorem [5j at any point x = 3>(y) of a suitable neighborhood where r ank(^£) = 
n — m, equation fp~T|) . by (fT6|) . is equivalent to 

FP^C/C/-^ = (J n _ OT , 0„_ m , m ) t^-ie = (19) 

from which we deduce that U' 1 ® = u «- m > m ^, where W is an arbitrary m x m polynomial 

V -'m / 

matrix, or B = U ( ® n ~ m ' m J py # The hyper-regularity of and the invertibility of \& immediately 
imply the hyper-regularity of W. But since W is square mxm and hyper-regular, it is unimodular. 



We introduce the notation: 



U = U[ ° n r m ' m ). (20) 



Lemma 2 For every Q G L — Smith ( Uj there exists Z G U m [-4] such that 



OB - ( Q Im )Z. (21.) 



J n—m,m 



Moreover, for every Q G L — Smith (Uj> the submatrix Q = (0 n _ m . m , I n -m) Q is equivalent to 
P(F), more precisely there exists a unimodular matrix L such that P(F) = LQ. 



Proof. The first part is an immediate consequence of Theorem [5] applied to U, which, by (|20p . is 

In 

Or, 



hyper-regular: there exist Q G ^-n[^] and R G U m [Jj] such that QUR = I m ). The result 

follows by setting R~ X W = Z G U m [|]. 

To prove the second part, we left multiply (|21 j) by (0 n _ m m , I n _ m ) and use ([18]) to get QU = 0. 
Then comparing with (|18p . where P(F)U = 0, and taking into account the hyper-regularity of [/, 
the existence of a unimodular matrix L such that P(-F) = LQ is proven. 



4.3 Integrability 



Let us denote by Qij = Ylk>o Qij^ * ne (*> i)"^ n polynomial entry of Q G L — Smith ^f/j obtained 
from Lemma [2j We also denote by to the m-dimensional vector 1-form defined by 



u)(x) 




{,1m,) m n — m ) Q{x)dx\ 



Xo 



J Xo 



\ 



(22) 
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— h (k) 

the restriction to Xq meaning that x G Xo satisfies L^F = for all k and that the dXj are such 

that dh\ x F = for all k in X . 

Since Q is unimodular, (I m , mjra _ m ) Q is hyper-regular of rank m. Thus, the forms u\, . . . , u; m 
are independent. 

Let us also recall that, if t\, . . . ,r r are given independent 1-forms in A 1 (Xo), the ^[^]-ideal T 
generated by n, . . . , r r , for an arbitrary integer r, is the set of all combinations with coefficients in 
of forms ?] A Tj with r/ an arbitrary form on Xo of arbitrary degree and i = 1, . . . , r. 

Note that if another set of independent 1-forms ki,...,k s of X is a generator of T, then s = r 
and there exists a unimodular matrix H G 1Xr[^j] such that r = .ff/c, where r = (ti, . . . ,T r ) T (the 
superscript T means transposition) and k = («i, . . . , K r ) T . 

Definition 7 Tl^e say i/iai £/ie ^[^]-idea/ T generated by T\, . . . ,r r is strongly closed if and only 
of there exists a matrix M G 1Xm[^] such that d(Mr) = 0. 

This definition is indeed independent of the choice of generators since if k is another vector of 
" ! ~ - — 11 f ™ H G U r [^|], and since d(Mr) = 0, we have d(MHn) = = 



d(M'K), with M' = MH G U r [^], which proves our assertion. 

Theorem 2 A necessary and sufficient condition for system (CD) to be flat at (xo,y ) is that there 
exist U G R — Smith (P(F)) and Q G L - Smith (tj^j , with U given by fity). such that the R[-j%]-ideal 
£1 generated by the 1-forms wi, ... , uj m defined by $22\) is strongly closed in Xq. 

Proof. Necessity: If system ([2]) is flat at (xo,y ), there exists $ = ((fo, (pi, . . .) meromorphic from 
a neighborhood of y i n ^oo to a neighborhood Xo of xq in Xo and bijective, such that x = &(y) 

o) 



implies F(ipo(y), (fi(y)) = and L\ : (F o $) = for all k > 0. According to Theorem [Q P(cp 



satisfies (fl~3j) . or (fT4"l) . and is hyper-regular: proceeding as in Proposition El its kernel P(y>o) _1 ({0}) 
cannot contain torsion elements since otherwise it would contradict the flatness property of the 
tangent trivial system (M™,0). Thus, in virtue of Lemma [U there exists W G U m [J^] such that 

pfao,^ = % ( \2 ,m ) *y ( 23 ) 

By Lemma [21 there exist Q G L — Smith ^{7^ and Z G U m [j] such that 

Ql_P(^o)| T ^ = ( } Z\- (24) 

Thus, left multiplying by (I m ,0 n _ mim ) and denoting by Q = (-^m, TO| „_ m ) Q, we get 

V (vo W) = i V (25) 



Taking the exterior derivative of x = fo(y) yields dx = P(ipo)dy, and thus, according to ([25]) . 

Z { X 1 Q\Jx = dy. (26) 
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The forms uj±, . . . , Lu m are independent, generate the ideal $7, and setting M = (Z)~ in (|26|) we get 
Muj = dy. Thus, taking the exterior derivative of both sides yields d(Moj) = 0, which proves that 
Q is strongly closed. 

Sufficiency: Assume that there exist U G R — Smith (P(F)) and Q G L — Smith ^{/^ such that 
Cl is strongly closed in a neighborhood Xo of xo in Xq. Let M € l-WJj] be such that d(Mui) = 0. 
Setting 7] = Muj, the 1-forms r]i,...,rj m generate f2, are independent in Xo since Q G L — Smith f f7 ^ , 

and, when expressed in the basis dx^ , depend on a finite number of terms whose coefficients depend 
on a finite number of derivatives of x. In the corresponding finite dimensional manifold, we have 
dr)i = 0, i = 1, . . . ,m, and, by Poincare's Lemma, there locally exists a mapping tpQ G C°°(Xo; M. m ) 
such that dipo = r/ = Muj = MQdx. In addition vfjQ is a meromorphic function of its arguments 
since its differential is, according to the previous relation. 

Denoting by y = ipo(x) for all x G X and * = (Vo, ;J^o, J^O, • • •) = (V'o, ipi, ip2, • • •), we have 
to prove that \P is a trivialization. 

Since Q G L — Smith (t/V writing as before Q = ( % ) with Q = (l m ,0m,n-m) Q and Q = 



Q 

(0n-m,m,In-m) Q, there exists R G Um[^] such that Qi7-R = I m and QUR = Q n - m . m . Setting 
W = RM~ l , we get QUW = M~ l and QUW = 0„_ mim , which, combined with dip = MQdx, 
yields QUWdifo = Qdx and QUWdipo = n —m- It results that there exists an element dz G ker(Q) 
such that UWdipo = dx + dz on Xq. But the latter relation, with the fact that dx is a basis of T*Xq 
and that, by the second part of LemmaEl we have Qdx = L _1 P(F)dx = 0, implies Qdz = 0, which 
means that dz G ker(Q) nker(Q) = {0}. We have thus proved that difto satisfies dx = UWdipo with 
W = RM^ 1 . If we denote by o"j the highest polynomial degree of the entries of the ith column of 
UW, we must have n < m + a\ + . . . + a m since otherwise this would contradict the surjectivity 
of UW, considered as the matrix H whose entries are the (J7W)^'s, mapping an open subset of 
+1 x ... x M°" m+1 to an open subset of M n . In addition, if we note a = max(<Tj|z = 1, . . . , m) 

dx 
9y° 



and y a = (y^ , ■ ■ ■ , y[ ai \ ■ ■ ■ , Vm , • • • , 2/m ) > S = J^? has rank n. Then the implicit system 



y = ipi {x) 



(27) 



has rank n with respect to x since its Jacobian matrix is a pseudo-inverse of S. Hence, by the 
implicit function Theorem, a local solution to (|27p is given by x = ipo(y, ■ ■ ■ , y^ a \ x, . . . , x^) for a 
suitable p G N. But differentiating tpg, using the fact that dF = 0, or equivalently P{F)dx = 0, 
and comparing with (fl8j) . we find that 990 is independent of (x, . . . ,x^), or x = <po(y). It results 
that <5 = (990 1 JjV^O) • • •) is the inverse trivialization of Vl/ which achieves to prove the Theorem. ■ 



Remark 2 T/ie strong closedness condition may be expressed in terms of distributions of vector 
fields by introducing the distribution D = {/ G TX|/J dF = 0, /J f2 = 0}, the symbol J denoting 
the interior product between vector fields and forms. This may be the subject of future works. 
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Before giving a characterization of strongly closed ideals, we need to introduce some notations 
and tools. First, let us denote by A P (X) the module of all the p-forms on X, by (A p (X)) m the 
space of all the m-dimensional vector p-forms on X, by (A(X)) m the space of all the m-dimensional 
vector forms of arbitrary degree on X, and by C q ((A(X)) m ) = C ((A p (X)) m , (A p+q (X)) m ,p > l), 
the space of all linear operators from (A p (X)) m to (A p+q (X)) for all p > 1, where C (V, Q) denotes 
the set of linear mappings from a given space V to a given space Q. 

We define the operator t), playing the role of the exterior derivative operator for unimodular 
matrices, by: 

d (H) k = cI(Hk) - HdK (28) 

for all m-dimensional vector p-form k in (A p (X)) m and all p > 1. Note that (|28p uniquely defines 
d (H) as an element of A ((A(X)) m ). 

We can prolong U for all /i G £ 9 ((A(X)) m ) and for all k G {A p (X)) m and all p > 1 by the 
formula: 

d(p) K = d(fi k) - (-l)V ck. (29) 
The oerator enjoys the following properties: 

Proposition 4 For a// /U G C q ((A(X)) m ), oil g £ N and a// k G (A p (X)) m ; w/zi/i p > 1 arbitrary, 
we have 

d(d(fi))n = 0. (30) 
In other words d 2 = 0, i.e. d is a complex. 

Proof. According to ([29]) . replacing p by t) (^) G ((A(X)) m ), we get 

t) (5 (^)) K = d (D (ft) k) - (-l) q+1 d (p) dK. (31) 
Since, again with (|29l) . d (D (p) k) = d 2 (p k) — (—l) q d (p dn) and, since d 2 = 0, we have 

d (0 (ji)) K = -(-l) q (t) (p) dK + (-l)V C?«) - (-l) q+1 d (p) dK 

= -{-l) q d(p) dn + (-l) q X) (p) dK = 
the result is proven. ■ 

We also have: 

Proposition 5 For all H G U m [j] and all p G C q ((A(X)) m ), with arbitrary q > 0, we have 

d{H)p + Hd (p) = d (Hp) . (32) 
In particular, if p = —H^^^d (H) G £\ ((A(X)) m ), we have 

(p) = p 2 . (33) 
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Proof. Let H G U m [±], k G (A^(3e)) m and fi G C q ((A(£)) m ), with q > arbitrary. ([28]) and 
([29D, yield 

T) (H) fin = d{Hfin) — Hd(fiK), d(fin) = d (fi) k + (— \) q fidn 

or 

5 (IT) UK = d{Hfi K ) - Hd (p) k - {-l) q Hfidn. 

In other words: 

(D (iT) fi + Hd{fi))K = d{Hfin) - (-l) q Hfidn =: (i2» k. 

This relation being valid for all k G (A^(X)| m and all p > 1, we immediately deduce ([32]) . 

If now /i = —H^ 1 ^ (H) G £i ((A^(3E)) ), we get —Hfi = d (H) and thus, according to what 
precedes, 

(IT/i) = d(H)fi + m (fi) = -d 2 (H) = 
or Hd (fj,) = —d (H) fi, or also t) (fjt) = —H^ 1 ^ (H) fi = fj, 2 , which achieves the proof. ■ 

A characterization of the strong closedness condition on Q is given by the next 

Theorem 3 The &[-rz]-ideal £1 generated by the 1- forms uj±, . . . ,uj m defined by $22\) is strongly 
closed in Xq (or, equivalently, the system (%,t%,F) is flat) if and only if there exists fi G 
C\ ((A(3£)) m ), and a matrix M G U TO [^] such that 

duj = fiuj, d(fi)=fi 2 , 5(M) = -M/i. (34) 

where we have denoted by uj the vector of 1- forms (uj\, . . . ,uj m ) T . 

In addition, if flff^[ ) holds true, a flat output y is obtained by integration of dy = Muj. 

Proof. If d(Muj) = 0, according to ([25]> . we have t) (M) u = -Mduj or duo = —M^d (M) uj. 
Setting fi = -M~ l D (M) G C\ ((A.^(je)) m ), which is equivalent to the last relation of ((HD, we have 
duj = —M~ 1 d (M) uj = fi uj. The conditions ([M]) are thus immediately deduced from ([33]) . 

Conversely, from the last identity of (I34p . we get — M^ 1 !) (M) = fi. Its combination with the 
first one yields duj = — M _1 t) (M) uj, or Mduj = —d (M) uj and, according to (I28p . we immediately 
get that d(Muj) = 0, i.e. Q is strongly closed. 

Finally, if there exists a matrix M G U m [J^] such that d(Muj) = 0, By Poincare's Lemma, there 
exist m functions yi, . . . , y m such that dy = Muj, which achieves the proof. ■ 



Remark 3 Condition {3$ may be seen as a generalization in the framework of manifolds of jets 
of infinite order of the well-known moving frame structure equations (see e.g. f^). 

From the necessary and sufficient conditions (134j) . an algorithm to test if a system is flat may 
be deduced. We proceed as follows. 
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Algorithm. 

1. We first compute a vector 1-form uj defined by (|22p . 

2. We compute the operator /i such that dw = fiu by componentwise identification. If no such 
At exist, the system is non flat. 

3. Otherwise, among the possible /u's, only those satisfying t) (//) = /x 2 are kept. If no fi satisfy 
this relation, the system is non flat. 

4. We then compute M such that d (M) = —M/j,, still by componentwise identification. 

5. Finally, only those matrices M which are unimodular are kept. If there are no such M, 
the system is non flat. In the opposite case, a flat output is obtained by the integration of 
dy = Muj which is possible since d(Muj) = 0. 

Remark 4 A non flat system is characterized by d(MuS) ^ for all M 6 U m [^-]. In this case, 
there exists a non zero m- dimensional vector 2-form t(M) such that d(Muj) = t(M) for all M G 
lt m [^]. Setting w = M~ 1 t(M), the 2-form w is clearly uniquely defined modulo (0,,dQ): indeed, 
if Mi and M2 are in U m [^] and if we denote by n = d(MiUj), W{ = M~ Ti and \ii = —Mf 1 d (Mi), 
i = 1,2, we have d(MiUi) = d (Mi) ui + Midui = T{,% = 1, 2, or doj = M^ 1 (r, — d (Mi) uj) = Wi + ^iUi, 
so that w\ — W2 = (fJ-2 — Ml) or w\ — vo<i = modulo (0, d£l), since, from A28\), it is easily seen 
that (a^2 — A^i) w i s a combination of elements of £1 and dVi. 
We thus have : 

Theorem 4 A system is non flat if and only if for every matrix M E ^m[^] there exists a non 
zero w(M) € (A 2 (X)) m such that: 

duo = n{M)u + w(M), dzu(M) = fi(M)w(M), d (fi(M)) = fi 2 (M) (35) 

with ijl(M) = -M^d (M). 

Proof. Since D (M) uj + Mdw = t(M) / for all M, we get du = M~ 1 t(M) + u{M)u) = 
w(M) + u(M)u after having noted M~ 1 t(M) = w(M), the first relation of 1135]) follows. By 
exterior differentiation, we get: = dzu(M) + t> (fi(M))uj — /i(M)du>. According to Proposition^ 
we deduce that t) (a(M)) = ^ 2 (M) (last relation of £35\) ) since, by definition, /x(M) = — M _1 c) (M). 
Thus = dw(M)+fj, 2 (M)-fj,(M) (zu(M) + h(M)uj), which proves the second relation of L35\). and 
(35)) is proven. 

Conversely, if \35\) is valid with w(M) 7^ for all M , we have d(Mui) = t) (M) uj + Mduj or 
d(Muj) = -Mn(M)u + M (zu(M) + u(M)u) = Mw(M) = t(M) ^ 
for all M € lt m [^:], which achieves the proof. ■ 

Therefore, one may introduce the equivalence relation "two systems are equivalent if and only 
if they have the same ideal f2 and their 2-forms w are equal modulo (O, dfi)", which suggests that 
it is possible to classify the non flat systems. In particular, a system has a non zero defect (see 
U6i[ ) if and only if it admits a non zero 2-form w modulo (£l,dQ), that might be interpreted as a 
generalized curvature (as opposed to flatness). 
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4.4 Some easy consequences 

We now show how several classical results of static feedback linearization |22l |2"U] , or in the case 
m = 1 |3 EH H2J H3] can be recovered as consequences of Theorem [2j 

If f2 is strongly closed, let us define o~% as the maximum degree of the entries of the ith column 
of P((fo), i = 1, ... ,m, which, according to Theorem^ locally yields: 



Definition 8 We say that a flat output y is minimal if N = YaLx a i * s minimal over all possible 
choices of U , Q and Z. 

Obviously, a minimal N always exists for flat systems. 

Corollary 1 A necessary and sufficient condition for system (Qjp to be static feedback linearizable 
(see [22, 20]) is that the strong closednes condition of Theorem^ holds true and that n = m + N 
for a minimal N . 

Proof. It is easily seen that ([2]) is static feedback linearizable if and only if ([2]) is L-B equivalent 
to a trivial system, the trivialization being such that </?q is a local diffeomorphism, which means 
that x is diffeomorphic to (y\, . . . , y( ai \ ■ ■ ■ , y m , ■ ■ ■ , y(°" m )) for a minimal y, i.e. n = m + N . ■ 

Remark 5 This Corollary generalizes the results of Jakubczyk and Respondek jl22f and Hunt, Su 
and Meyer ]2(tf in a twofold manner: first, it is not restricted to affine systems and second it 
applies not only in a neighborhood of an equilibrium point but of any trajectory around which the 
variational system is controllable. 

Corollary 2 If m = 1, a necessary and sufficient condition for flatness is that £1 is closed in the 
ordinary sense, i.e. duo\ = r Awi for some 1-form t. Furthermore, the system is flat if and only if 
it is static feedback linearizable. 

Proof. If 17 is generated by a single 1-form u±, the unimodular matrix M must be a non zero 
element of 8. and the strong closedness of O reduces to ordinary closedness. Thus the system is 
flat if and only if Q is closed. Assuming closedness, by Frobenius' Theorem, we deduce that there 
exists a scalar function ipo such that dipo = Mux an d that (|36p holds true with N = n — 1. Using 
Corollary [21 the system is static feedback linearizable. The converse is trivial since every static 
feedback linearizable system is flat. ■ 



5 Examples 

5.1 Non holonomic car 

Consider the 3 dimensional system in the x — y plane, representing a vehicle of length I, whose 
orientation is given by the angle 0, the coordinates (x, y) standing for the position of the middle 




(36) 
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of the rear axle, and controlled by the velocity modulus u and the angular position of the front 
wheels ip. 

x = u cos 9 

y = u sin 9 (37) 
9 = j tan ip 

Since n = 3 and m = 2, n — m = 1 and (|37|) is equivalent to the single implicit equation obtained 
by eliminating the inputs u and (p: 

F(x, y, 9, x, y,9) = x sin 9 — y cos 9 = (38) 

We immediately have: 

P(F) -(— + —— — + —— — + — — 

\dx dx dV dy dy dt' 89 dt ) dt / ,., <v 

d n d N loJ) 

sin^— , — cos# — , x cos u + y sin I 
dt dt y 

Setting E = xcos9 + ysm9, we apply the Smith decomposition algorithm of the Appendix: moving 
the last column (of degree zero) to the first place by a permutation with the two others of degree 
one in any order, we see that rank (P(F)) = 1, that P(F) is hyper-regular and we get 

/ 1 
U= 1 

\ 1_ cos 9 d sin 9 d 

\ E E dt E dt 

Thus 

u = u( 0l > 2 ) = ( 1 

\ h J \ cose d sing d 
\ E dt E dt 

with I2 the identity matrix of M 2 . Again, computing Q G L — Smith ( U ) yields 

1 
Q = I 1 

sine d cos 9 d -1 

E dt E dt 

I dx\ 

Multiplying Q by the vector \ dy J , the last line reads 

\d9 J 

-g (sinOdx — cos 9dy + (a: cos 9 + y sm9)d6) = ^d{xsm.9 — ycos9) and, by (f38]) . identically 
vanishes on Xq. 

The remaining part of the system, namely J ||^^(iy^=^ a;i ^is trivially strongly 

closed with M = I2, which finally gives the flat output y\ = y and y2 = x. We have thus recovered 
the flat output originally obtained in |41L I40j . up to a permutation of the components of y. 
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5.2 Non holonomic car (continued) 



Other decompositions of P(F), given by (|39|) . may indeed be obtained, but they are all equivalent 
in the sense that one decomposition may be deduced from another one by multiplication by a 
unimodular matrix. However, the resulting vector 1-form uj, contrarily to what happens in the 
previous example, may not be integrable. Our aim is here to show how the generalized moving frame 
structure equations may be used to obtain an integrable Mu. Such an example is provided by 

/ cos6> 

restarting the right-Smith decomposition of P(F) by right-mutiplying it by I sin# 1 | and 

\ 1 

using the formula sin#^(cos#) — cos#^(sin#) = —6, we obtain 



U 



' COS I 

sin/ 

V 



2 cos 2 04 
| sin 9 cos 6 




dt 



i (x cos 9 + y sin 9) cos ( 
-* (x cos 9 + y sin 9) sin ( 



U 



( 



\ 



cos 



I (x cos 9 + y sin 9) cos 9 
1 — i sin 9 cos 94t j{x cos 9 + y sin 0) sin 



and then QgL - Smith (u): 



Q 



— tan 9 


i sin0 cos #4 



1 



W0| 





1 



OA -i(xcos9 + ys'm9)cos9 



The vector 1-form u> = (uji,u>2) t is obtained by multiplying the two first lines of Q by 
(dx , dy , d9) T : oj\ = —tan9dx + dy and 002 = d6. We have duj\ = — — g g d0 A and cL>2 = 0. 





Hence, in order that dw = /j,uj, a possible choice is \x - 
arbitrary meromorphic function. One can verify that /j, is admissible if rj = 
t) (n) = = 11 2 . Again, by componentwise identification, one finds M = 
have 

( —tan9dx %sd0 + dy \ ... r . 



where 77 is an 



-2 ZS1 3 ^ to ensure that 

cos y 

1 



cos^ 

1 



and we 



Thus, setting 



dyi 
dy2 



Muj, 



yi 



dd 



y — x tan 9, y2 



which is another possible flat output: it is easily checked that the inverse L-B isomorphism is given 
by x = -|i cos 2 y 2 , y = y\ - siny 2 cosy 2 , = y 2 . 
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5.3 The pendulum 



We consider a pendulum in the vertical plane studied in [IB], of length I and inertia J, whose mass 
m is concentrated at its end point C. It is controlled by the two components of the force F applied 
to the opposite end point A of the pendulum. Introducing an inertial frame (0, x, z), it is modelled 
by 

X = U\ 

z = u 2 (40) 
aO = —u\ cos 9 + (ii2 + 1) sin 9 



with x = z = ?£L, (xc, zc) being the coordinates of C, 9 being the angle between the pendulum 



XC ~ — *C 
9 ' 9 '' 

and the vertical axis, u% = u?, = ^ — 1, (F X ,F Z ) being the components of the force F, and 

a = mg~r 

An implicit model is given by: 



F( 



a9 + x cos 9 — (z + 1) sin I 



0. 



(41) 



Though this is a second order system, it can be easily transformed into a first order one by setting 
x = v x , z = v z and 9 = vg. We thus obtain the 4 dimensional implicit system 



x — v x = 

z — V z = 
9 - Vg = 

at># + v x cos ( 



(42) 



(v z + 1) sin 6» = 0. 



However, it can be easily verified that all the results of this paper can be extended word for word 
to higher order systems. Thus, because of its smaller dimension, we prefer using (|41h instead of 

T2D- 

The variational system corresponding to (|4ip is given by 



cos 



d 2 
)— 

dt 2 ' 



d 1 



SmU df^ 




(43) 



where b = 'x sin 9 + (z + 1) cos 9 
Using the identity ^cos6>^ 
setting E = a9 — b, we have 



-a cos I 



sin 



) d2 
dP 



a sin ( 



+ 



-4 

'dP 



6 (1) 



b and 



cos 



sm 




1 sin 9 d 2 


cos 9 


E E It? 


E 


1 








1 



d 2 



dt 2 



(1,0,0) 



which proves the hyper-regularity of P(F) and we get 



U 



a sin 9 cos 8 d 2 a cos 2 e d 2 



E 

a sin 9 



E 



9 E 9 ^ 

a sin^ 9 d z L -i 

"~ E~~ ST," 1 " 1 

sing ar 
E dF 



a sin 6* cos d 
E IP 
cos 6 ar 
-E dl 2 " 



C7 



a sin 9 cos d 2 
E dF 



sm & 
\ E 



E dP n 
sine d 2 
E dP 



+ 1 



d 2 



+ 1 



E dP ., 
a sin 8 cos d 



E 
cos 6 d 
£ dP 



,//- 
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Now left decomposing U, we get 



Hence, 



and 



Q 



dx 
dz 
d6 



1 















1 





sin 9 


d 2 


cos 9 d 2 


1 


E 


w 


E df 2 



Q 



( dz — a sin 6d6 

dx + a cos 9d6 

y -gj ^cos 6dx — sin Od'i + a<i(9 






1 


— asin^ \ 








1 





acos# 


IN 













1 J 




V o 


o / 






1 




—as'mO 


1 







a cos 


cos 9 d 2 


sin 


d 2 


ad 2 b 


E dt z 


E 


IF 1 


EdP E 



bd6 



dz — a sin 6d9 
dx + a cos ( 




Thus, the strong closedness condition holds true: setting M = I2, we obtain 

d(z + a cos 9) = dyi, d(x + a sin 6) = di/2 



d(z + acos#) 
d(x + asin^) 




or 

y% = z + a cos 0, 7/2 = ^ + sin 

which represents, up to a permutation of yi and 1/2, the coordinates of the Huygens oscillation 
center, already found in [26l fl6] . 



6 Concluding remarks 

We have proved that flatness is equivalent to the strong closedness of the ideal of 1-forms rep- 
resenting the differentials of all possible trivializations, using an approach which is invariant by 
endogeneous dynamic feedback. Moreover, we have separated the algebraic characterization of 
their differentials (Lemmas [T] and [2]) and the integrability aspects, that may be seen as a gener- 
alization, in the framework of manifolds of jets of infinite order, of the well-known moving frame 
structure equations. The computation of flat outputs in non trivial examples show the applicability 
of our results. Note that the solutions to these examples were already known since long. However, 
some other classes of examples, whose solutions are not presently known, are in preparation. 

We have chosen to present our results in terms of differential forms, in contrast with a large 
part of the control literature in this domain, where vector fields and Lie brackets are prefered. 
The main argument in favor of the former language is that the results are easy to state and are 
constructive: flat outputs are directly computed by integration of the strongly closed ideal £1. It 
might be possible however to express them in terms of distributions of vector fields and this could 
be the subject of future works. 

Finally, our necessary and sufficient conditions might open a way to a classification of non flat 
systems by the introduction of a 2-form that may be interpreted as a generalized curvature (see 
Remark |4|) . 
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Appendix 

A The Smith decomposition algorithm 

We consider matrices of size p x q, for arbitrary integers p and q, over the principal ideal do- 
main &[■%], here the non commutative ring of polynomials of ^ with coefficients in the field & of 
meromorphic functions on a suitable time interval 3. The set of all such matrices is denoted by 
For arbitrary p 6 N, the set of unimodular matrices of size p x p is the subgroup 

of M PiP [^] of invertible elements, namely the set of invertible polynomial matrices whose inverse is 
also polynomial. 

The following fundamental result on the transformation of a polynomial matrix over a principal 
ideal domain to its Smith form (or diagonal reduction) may be found in \10\ Chap.8]): 

Theorem 5 Given a (fixv) polynomial matrix A over the non commutative ring with [i<v, 

there exist matrices V G and U E U^^] such that VAU = (A,0) where A is a /i x /i (resp. 

v xv) diagonal matrix whose diagonal elements, (Si, . . . , 8 a , 0, . . . , 0), are such that 5i is a non zero 
■^-polynomial for i = 1, . . . , a, and is a divisor of Sj for all a > j > i. 

The group of unimodular matrices admits a finite set of generators corresponding to the follow- 
ing elementary right and left actions: 

• right actions consist of permuting two columns, right multiplying a column by a non zero 
function of .ft, or adding the jth column right multiplied by an arbitrary polynomial to the 
ith column, for arbitrary i and j; 

• left actions consist, analogously, of permuting two lines, left multiplying a line by a non zero 
function of ^, or adding the jth line left multiplied by an arbitrary polynomial to the ith 
line, for arbitrary i and j. 

Every elementary action may be represented by an elementary unimodular matrix of the form 
Tij(p) = I v + lijp with ljj the matrix made of a single 1 at the intersection of line i and column j, 
1 < z , j < , and zeros elsewhere, with p an arbitrary polynomial, and with v = m for right actions 
and v = n for left actions. One can easily prove that: 

• right multiplication ATij (p) consists of adding the ith column of A right multiplied by p to 
the jth column of A, the remaining part of A remaining unchanged, 

• left multiplication Ti j j(p)A consists of adding the jth line of A left multiplied by p to the ith 
line of A, the remaining part of A remaining unchanged, 

• Trj(p)= Ti j(-p), 
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• Tij{l)Tj,i(- l)Tjj(l)A (resp. AT^l)!)^-!)!^!)) is the permutation matrix replacing the 
jth line of A by the ith one and replacing the j th one of A by the ith one multiplied by — 1 , 
all other lines remaining unchanged (resp. the permutation matrix replacing the ith column 
of A by the jth one multiplied by —1 and replacing the jth one by the ith one, all other 
columns remaining unchanged). 

Every unimodular matrix V (left) and U (right) may be obtained as a product of such elementary 
unimodular matrices, possibly with a diagonal matrix D(a) = diag{«i, . . . , a v } with «j £ ^, «j ^ 0, 
i = 1, . . . , v, at the end since Tij(p)D(a) = D(a)Tij(-^-paj). 

In addition, every unimodular matrix U is obtained by such a product: its decomposition yields 
VU = I with V finite product of the Tjj(p)'s and a diagonal matrix. Thus, since the inverse of any 
Tij(p) is of the same form, namely Tjj(— p), and since the inverse of a diagonal matrix is diagonal, 
it results that V^ 1 = U is a product of elementary matrices of the same form, which proves the 
assertion. 

The algorithm of decomposition of the matrix A consists first in permuting columns (resp. 
lines) to put the element of lowest degree in upper left position, denoted by a^i, or creating this 
element by euclidian division of two or more elements of the first line (resp. column) by suitable 
right actions (resp. left actions). Then right divide all the other elements a± jk (resp. left divide 
the a^i) of the new first line (resp. first column) by a\^. If one of the rests is non zero, say r\ jk 
(resp. rk,i), subtract the corresponding column (resp. line) to the first column (resp. line) right 
multiplied (resp. left) by the corresponding quotient q\ k defined by the right Euclidian division 
ai,k = ai,iQi,k + n,fc (resp. q kjl defined by a kA = q kjl a 1;1 + r k>1 ). Then right multiplying all the 
columns by the corresponding quotients q± >k , k = 2,...,v (resp. left multiplying lines by q kj i, 
k = 2,...,fi), we iterate this process with the transformed first line (resp. first column) until it 
becomes (fli i, 0, . . . , 0) (resp. (ai,i, 0, . . . , 0) T where T means transposition). We then apply the 
same algorithm to the second line starting from 02,2 and so on. To each transformation of lines and 
columns correspond a left or right elementary unimodular matrix and the unimodular matrix V 
(resp. U) is finally obtained as the product of all left (resp. right) elementary unimodular matrices 
so constructed. 
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